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Following a recent work by M. Hillery et al. [ Phys. Rev. A 68, 032314 (2003)] we introduce 
a scattering model of a quantum random walk (SQRW) on a hybercube. We show that this type 
of quantum random walk can be reduced to the quantum random walk on the line and we derive 
the corresponding hitting amplitudes. We investigate the scattering properties of the hypercube, 
connected to the semi-infinite tails. We prove that the SQRW is a generalized version of the coined 
quantum random walk. We show how to implement the SQRW efficiently using a quantum circuit 
with standard gates. We discuss one possible version of a quantum search algorithm using the 
SQRW. Finally we analyze symmetries that underlie the SQRW and may simplify its solution 
considerably. 
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I. INTRODUCTION 

Quantum random walk is a theoretical concept con- 
ceived to simulate certain algorithms using quantum me- 
chanical elements, i.e. unitary operators and measure- 
ments I n particular, it has been shown recently that 
it is possible to use a quantum random walk to perform a 
search in a database with the topology of the hypercube 
faster that it can be done classically It is an oracle 
based algorithm, which is optimal in its speed. Another 
successful application of quantum random walks has been 
demonstrated by Childs, et al. citechilds2003 who have 
also constructed an oracle problem that can be solved by 
a quantum algorithm exploiting quantum random walk 
exponentially faster than any classical algorithm. 

These two examples justify the general hope that quan- 
tum random walks might be able to solve some problems, 
based on random processes (e.g. Monte Carlo meth- 
ods, 2-SAT, graph connectivity, etc.) faster than corre- 
sponding classical algorithms. 

In general, there are essentially three types of quan- 
tum random walks. Firstly, let us mention the so-called 
coined quantum random walk (CQRW), that is a discrete 
time walk which makes use of an additional quantum sys- 
tem, the coin [I| . The second type of the quantum ran- 
dom walks is described by a continuous (Hamiltonian) 
dynamics of a quantum system |6j. The third type of 
quantum random walks based on physical model of opti- 
cal multiports has been recently been introduced by M. 
Hillery et al. 0,0. 

Before we proceed we note, that a quantum random 
walk as discussed by D. Aharonov et al. 11] is basically 
a mapping ip : Q\> — ► C d , which is updated at each step 
by a function tp(x) i— > F[ip(y) : (xy) e Qe], where Q — 
(Gv-i Gs) is a graph with vertices Gv and edges Gs- 

We therefore can say that quantum random walk is a 
special instance of the quantum cellular automaton [Tol 
Il3 ] . The classical cellular automaton is a concept general 
enough to accommodate virtually any algorithm; more 



precisely, any Turing machine can be simulated using a 
cellular automaton. 

The CQRW is usually denned on regular graphs (each 
vertex having the same number of outgoing edges). The 
definition on non-regular graphs is also possible, and 
some interesting algorithms are based on this version 
[T3 | . However, the latter version does not possess the 
symmetries of the former one, nor its neat tensor prod- 
uct structure (the unitary evolution operator CQRW on 
the regular Cayley graph commute with generators of the 
underlying group). Instead, the whole graph must be ad- 
dressed, by means of an oracle which tells us whether any 
two vertices are connected by an edge [l5|. which causes 
a considerable growth of the resources. 

In this paper we will focus our attention on a quantum- 
optical model of multiports 0, @ which describes a pos- 
sible physical realization of specific quantum random 
walks. In this scheme, we have an array of multiports 
(see e.g. [T^ ] and references therein), interconnected with 
optical paths. A photon is launched into one path and 
is transformed by the action of the multiports. This ac- 
tion can be described as a scattering process, therefore, 
we will refer to this scheme as the scattering quantum 
random walk (SQRW). 

The SQRW is more viable from the experimental point 
of view, can be extended to non-regular graphs and is 
equivalent to CQRW on the regular graphs. 

We will investigate a particular arrangement of the 
multiports when are localized at the vertices of a hyper- 
cube. The array of multiports effectively acts as a scat- 
tering potential, when connected to semi-infinite tails. It 
can be endowed with two characteristic values, the reflec- 
tion and transmission amplitude for photons. This was 
done in Ref. [3, |8| f° r a special two-dimensional hyper- 
cube. 

Our paper is organized as follows: in Sec. [H]we define 
the SQRW on the hypercube. In Sec. IIHI we show how 
the SQRW on the hypercube may be reduced to the cel- 
lular automaton on the line. In addition we will compute 
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the hitting amplitude and we will make some other sim- 
ulations. In Sec. II VI we will investigate scattering prop- 
erties of the hypercube connected to semi-infinite tails. 
Sec. I VIII contains the proof that the SQRW is equivalent 
to a generalized version of the coined quantum random 
walk. In Sec. IVII we show how to implement the SQRW 
efficiently on the quantum circuit. In Sec. IVIIII we dis- 
cuss one possible version of search algorithm using the 
SQRW. Finally in Sec. IIXI we will analyze symmetries 
that underlie the SQRW and which may considerably 
simplify the solution of the model. 



II. THE DEFINITION OF SQRW 

The SQRW was first presented in 0. The technique 
behind its implementation are the multiports ^(|: lin- 
ear optical elements, interconnected with optical paths. 
Each multiport has in general different number of inputs 
(which at the same time serve also as outputs). A co- 
herent superposition of photons entering the multiport is 
transformed into another coherent superposition of pho- 
tons outgoing from the multiport. The multiports are 
the vertices of a graph C/y, and the optical paths are its 
edges. 

The photon travelling between multiports x, y is de- 
noted \xy). The Hilbert space on which the state of the 
photon is defined is spanned by vectors |xy), {xy) £ Ge 
and can be decomposed into direct sum of Hilbert sub- 
spaces H = ® xeGv Ti x where H x = span{|ya;) : (yx) £ 
Gs}- That is, the base of Tt x are the states of photon in- 
going into multiport at x. For convenience we introduce 
the Hilbert subspace spanned by the states of photon out- 
going from the vertex at x: H x = spa,n{\xy) : (xy) £ Gs} 



The evolution operator of states in Ti is U = ® x 



u x 



where U x is isometry TL X — > TL X (onto, hence linear, hence 
unitary). Since the subspaces are orthogonal, U is uni- 
tary. 

The concrete realization of unitary operator U x reflects 
the fact that the multiport partially reflects and partially 
transmits the ingoing photon. Denoting the reflection 
and transmission coefficients respectively r and t we have 



U x \yx) = r\xy) + t \xz) 



(1) 



where (xy) £ Ge For U to be unitary, these coefficients 
must satisfy the relations 



\r\ 2 + (d-l)\t\ 2 = 1 
(d-2)\t\ 2 +r*t + rt* = 



(2) 
(3) 



The operator U x (for any x) has in the natural basis the 
matrix the form (U x )ij = rdij + t(l — <5,j). The rela- 
tions (QJ) may be satisfied by setting by u sing the Grover 
coefficients r = § — 1, t = |, see (17| . The pseudo- 
eigensystem of any such U x can be readily computed 



( "pseudo" means that we set an isomorphism between Ti x 
&ndTL x s.t. \yx) = \xy)). Since U x — 2\s)(s\ — 1 (|s) is the 
complete superposition over the basis of the domain of 
U x , which we denote |1), . . . , \d)), the eigenvectors are |s) 
(with pseudo-eigenvalue 1) and linearly independent (but 
not orthogonal) set {|1) — |fc) : k = 2, . . . , d} (with (d— 1)- 
degenerate pseudo-eigenvalue -1). Making the direct sum 
of these eigenvectors gives us a pseudo-eigensystem of U, 
with pseudo-eigenvalues ±1. 

There are many other choices of the reflection and the 
transmission coefficients, one set of them is the following 
(we will use it when necessary): 



1 

dP 



t 



r= \ 1- 



(4) 



d?P 



where cos 6 



If we set 



,]~ d/2 with p > 1/2. 
p = 1, then the reflection amplitude goes to lim^oo r = 
— i + z-^. Obviously, for large dimensions, almost all 
of photons will be reflected. The multiports with these 
coefficients will be called symmetric multiports. We will 
later compare both sets of coefficients with respect to 
their mixing properties. 

From now on we will be dealing with the multiports 
arranged into the form of the d dimensional hypercube. 
The edges of the hypercube will form two-way optical 
paths. 

The (i-dimensional hypercube is the Cayley graph G = 
(Zj, [d]) where [d] is the set of generators of the additive 
(mod 2) group 1^ (the binary strings with only one 1). 
For any x,y e Zj we set the scalar product xy — x\y\ + 
• • ■ + x c iyd (mod 2). The norm |a;| = ^pxx is the Hamming 
weight (the number of Is in x). The set £ w = {x : \x\ = 
w} is called the Jayer of the hypercube. Obviously [d] = 
ii. 

For simplicity, we will denote the basis states of H for 
the (i-hypercube \xy) like \x; a) where x is the vertex and 
a = 1, . . . , d is the generator s.t. x + a = y. 



III. THE TOPOLOGY OF THE HYPERCUBE 

The hypercube may be broken up in individual lay- 
ers t w , i.e. sets of vertices with equal Hamming weight. 
There is a special class of vectors from H, which is 
closed w.r.t U and whose members may be described by 
fewer coefficients, thus simplifying the evolution equa- 
tions. Namely, these are the vectors s.t. (x;a\ip) 
is the same for all \x\ = w and for all a. Under this as- 
sumption, the vectors are specified by coefficients {ip Wl ±} 
where ip w .± — (x; a\ip) with \x\ = w, \x + a\ = w ± 1. 

The reduced equations for evolution of the coefficients 
"0u>,± are given from the assumptions that each vertex 
from £ w has a fixed number of edges which connect it 
to the previous and next layer (with Hamming weight 
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FIG. 1: Action of the multiport (U) on the ingoing photon 
(in) on the dim 3 hypercube. Coherent superposition of 3 
photonic excitations is created. 



The classical and quantum hitting probabilities are re- 
lated like 




5 10 15 20 



having w ± 1). We note that a vertex \x\ = w is con- 
nected with w edges from previous layer and with d — w 
edges with next layer. Since the coefficients assigned to 
the projection to a given layer are fixed, we get the 
recursive relations: 



(Ulp) w<+ = tW1p w -i t+ + [t(d -W—l) + 7#tu+l, 

(Uip) Wl - = t(d ~ w)ip w+ i- + [t(w - 1) + r]^u<-i,+ 



(5) 



For w = 0, d these equations still hold wherever it makes 

sense, i.e. for ^o,-H tyd, The coefficients ip-i,., i^d+i,- 

are neglected as long as they are multiplied by zero in 
the equation. 

The evolution which is governed by these equations is 
called the symmetric SQRW. 

If the initial state is = then we immediately 
obtain the expression for the hitting amplitude 



lk,_(d) = [t(d-l) + r] (d-1)! 



d-l 



Vd 



(6) 



Using the Stirling approximation log n! ss n log n — n we 
can express for large d the above equation as ipd,-(d) = 
6(e~ d ). We see that the hitting amplitude drops expo- 
nentially with d. Classically, the probability that we 
are at a given vertex from £ Wl is p w , for which holds 
(Wp) w = ^wPw-i + \ {d-w)p w+ i - 1. From this we get 
the hitting probability (W d po) = -gzd\. By simple means 
we show that quantum and classical probabilities to get 
from the vertex Ire j = to vertex \x\ = d in d steps are 



(e/2) 



Vd 



(7) 



FIG. 2: The ratio of hitting probability for classical and quan- 
tum random walk on the hypercube, related to the dimension 
d. 

to simulate coined quantum random walk on the segment 
of line with coin dependent on the position. The state of 
the coin is described by vector [Vv-H VV-L w = 0,...,d 
and is updated at each step with update rules ©. We 
might be troubled by the fact that the update rule is not 
unitary. But we need the unitarity only to conserve the 
inner product {a\a). Actually, the inner product 



(v¥> = E 



w=0 



(|Vv + | 2 + IW-l 2 ) 



(8) 



is conserved. 

Though we have simplified the problem by the assump- 
tion of symmetric initial values, we are still far from its 
explicit solution. The solution would rely on the path in- 
tegration along different paths by which two sites can be 
connected in a presupposed number of steps. Each path 
would be assigned a complex amplitude (basically some 
product of r, t) , and by adding all the relevant paths to- 
gether, we would get the amplitude distribution over the 
hypercube. This normally gives us enormous combinato- 
rial expressions, which are difficult to interpret. 

Since the SQRW on the hypercube with symmetric ini- 
tial states is equivalent with one-dimensional (quantum) 
random walk on the finite sequence of layers £ w , we may 
explore the probability distribution p n (w) over the layers 
for an initial state |V>o) = Yla "73^ ' ' ' ^' a )' wnere Pn{w) 
is the expectation value of the operator 



M m := 



E \x;a)(x;a\ 



(9) 
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i.e. p n (w) = (i/} \(Ut) n M w U n \ijjo). We simulate the evo- 
lution of p n (w)\ see Fig|3 

p(w) 



. 8 
. 6 




FIG. 3: The probability distribution of the SQRW, p n {w), 
for the hypercube of the dimension d = 50, with a sym- 
metric initial state localized at the vertex ... 0, and steps 
n = 1, . . . , 100. The multiports are Grover. 

We see that the SQRW on the hypercube is effec- 
tively isomorphic to a dynamics of a resonator: we start 
with one excited site, and the excitation propagates as a 
Gaussian packet along the resonator. As the packet hits 
the boundary, the corresponding site is excited, and the 
packet is reflected in the opposite direction. 
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FIG. 5: The probility distribution of the symmetric SQRW on 
the hypercube of dim 50, for steps n to 250, with the initial 

state tpd/2+i,± = Vy^Ci/Jfi) (other tp-s=0) and symmetric 
multiports. 




w 



FIG. 6: The probability distribution of the symmetric SQRW 
on the hypercube of dim 50, for steps n to 250, with the 
initial state ipo,+ = V'd,- — (other tp-s=0) and Grover 
multiports). 



FIG. 4: The probability distribution of the symmetric SQRW 
on the hypercube of dim 50, for steps n to 250, with the 
initial state ipo,+ = ipd,- = -^=g (other ip-s=0) and symmetric 
multiports. 

From Fig. [3] we see that the hypercube acts as a res- 
onator when it comes to the evolution of the probability 
distribution over the layers. 

We have performed detailed simulations for the sym- 
metric SQRW with initial conditions such that either- 
only the coefficients ^0,+ ; ipd.- are non-zero, or ipd/2,± 7^ 
0. There are two choices of the multiports: the Grover 
and symmetric multiports. We make the simulations for 



both of them, and for both initial conditions (see Figs. |H- 
EJ. We see the periodicity of the evolution: another fea- 
ture akin to the resonant behavior. Note that the Grover 
multiport has much better mixing properties than the 
symmetric multiport, owing to the fact that it is more 
distant from the unity. 



IV. HYPERCUBE AS A SCATTERING 
POTENTIAL 

In Ref. a model of 2-d hypercube with semi-infinite 
tails attached to the vertices with Hamming weights 
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FIG. 7: The probability distribution of the symmetric SQRW 
on the hypercube of dim 50, for steps n to 250, with the ini- 
tial state y>d/2+i,± = Vy2(d/^fi) (other ^-s=0) and Grover 
multiports. 



and 2 has been studied. Each tail has been supposed 
to be a 1-d lattice with perfectly transmitting multi- 
ports. Along one tail, a photon enters the hypercube, 
and emerges on the other side. It is possible to calcu- 
late explicitly the transmission coefficient of the whole 
structure. In the present section we will study analogous 
problem for an arbitrary-dimensional hypercube. We will 
utilize some symmetry assumptions, which let us perform 
the calculation (or at least the simulation) for arbitrary 
higher dimensions. 

The scheme we consider is shown on the figure [SJ 
Now everything is as before, except that the multiports 
at the vertices with Hamming weight and d have reflec- 
tion/transmission coefficients r — ^xj — 1 1 = ^tj- The 
multiports outside the hypercube are perfectly transmit- 
ting. The initial state is a photon travelling between 
vertex —1 to vertex |x| = from the hypercube. The 
state of the whole system is described by complex num- 
bers a w> ± which represent the amplitude that the photon 
travels from a vertex \x\ = w and is directed to the next 
or previous layer l w ±\. Now we let w running through 
Z. The resulting relations are: 

(Ety) ,+ = + [(d - 1)* + f]ipi,- 

{Uip) - = r0_i )+ + dtipi _ 

(Uip) d ,+ = dtipd-i,+ 

(Uip) Wt+ = twip w _ li+ + [r + (d-w- l)t]ip w+1 - 

(U i>) w - = t(d - w)ip w+1 - + [t(w - 1) + r]^ w -i,+ 

(10) 

We begin with the particle in the state | — 1, 0), i.e. a 



FIG. 8: The scattering potential (hypercube dim 3). The 
vertices outside the hypercube are denoted — oo, . . . , —2 and 
d + 1, . . . for hypercube of dimension d. 



particle localized at the vertex just left of the hypercube 
on the tail, and pointing to the right (see Fig[HJ|. 
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FIG. 9: The scattering probability of dim 10 hypercube, for 
a photon incoming from the source (S) (n is the number of 
steps). The multiports are symmetric. 



We have simulated the probability that a particle in- 
coming from the left will be absorbed by the detector 
after n steps (see Fig. |5J). This means that the system 
evolves for n steps, and then a projection on the basis 
{\d, +)} is made. The result shows periodic beats of the 
probability of absorption of the photon by the detector. 
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V. NON-SYMMETRIC INITIAL STATE 

When wc impose a symmetry condition on initial states 
the whole problem becomes linear in the dimension of the 
hypercube (normally its complexity is exponential in d). 
But we also may be interested in the behavior appears 
as a result of phase differences in the components of the 
initial state. Given the d-hypercube with semi-infinite 
tails attached to the vertices ... and 1 ... 1, and the 
initial state 



l^>=I>j|0...0;i) 



(11) 



Now the amplitude of U d \tfjo) to project on |1 . . . 1; +) is 
given by the sum of amplitudes to traverse from ... 
to 1 ... 1 in d steps, which is t^^t. For initial state to 
be a concrete |0 . . . 0; j) we have (d— 1)! such paths. The 
situation is analogous for all j , with each initial direction 
j contributing the factor to the amplitude. The overall 
amplitude is 



d 

(l...l,+|L/<Vo>=£7^- 



l)!^" 1 ^ (12) 



The probability of detecting the particle at the state 
|1 . . . 1; +) depends only on Y^j=i 7r l n this sense the hy- 
percube with tails attached behaves like Mach-Zchnder 
interferometer. 



VI. IMPLEMENTING THE SQRW 

Until now we haven't discussed the question whether 
it is feasible to implement the SQRW. To build a whole 
network of multiports we need exponentially growing re- 
sources (the number of vertices grows exponentially). 
However, to encode the states we need only d|~logef| 
qubits. So we can ask a question: Is it possible to build a 
network of quantum gates operating on the qubit register 
of this size? This is most easily done only on the hyper- 
cube without the semi-infinite tails attached; however, 
it is also possible to implement this scheme, by adding 
some overhead of gates to the network. We need d qubits 
for the position register |x) and at least flog rf] qubits for 
the direction register \<f>). The first part of one applica- 
tion of the unitary operator U is controlled negation of 
each bit of x depending on \ip). The second part is the 
transformation of the state \<p), so that the action of the 
multiports is accounted. More precisely, the first part is 



and the second part reads 



|s>|o>->|s>[r|a>+X>>], 



(13) 



(14) 



for each a = 1, . . . , d. The first part described by Eq. (|13|) 
can be implemented using variant of CNOT gate. The 
CNOT gate operates on two qubits such that it negates 
the first (target) qubit, iff the second (control) qubit is 
nonzero. In the symbolic form, we have: CNOT = a x <£> 
|1)(1| + 1 ® |0)(0|. We employ the 0CNOT gate, which 
differs from the CNOT gate in that it has a ci-dimensional 
control state, unlike a single qubit. If the control state is 
\<j>) (the accepting control state), then the target qubit is 
negated, otherwise, it is kept in the original state. The 
operational form of the </>CNOT gate is 



(/.CNOT = a x 



1®(1-|0)<0|) 



(15) 



Obviously, the 0CNOT is unitary. 

The operation l|13fl can be implemented using d 
0CNOT gates (see Fig lTUfl . Each gate operates on a dif- 
ferent qubit from the position register. If the gate op- 
erates on the a-th qubit, the accepting control state is 
chosen to be \a) (see Fig. ITU)) . 
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FIG. 10: The gate which implements the SQRW on the 
d-hypercube. The input state is the position register (d 
qubits labeled 1, . . . , d) and the direction register \ip). There 
are d </)CNOT gates stacked together, with accepting states 
|1), . . . , \d), which change the position register, and the gate 
M which changes the direction register. 

The operation i|14|) is implemented using a single uni- 
tary operation M operating on the direction register. It 
corresponds to the transformation of the state due to the 
multiports. Its reads 



M = £[r|a)<a|+i£ 
The d x d matrix form of M is 



(16) 



M 



(r t 
t r 



V 



t r) 



(17) 



Consequently, the unitary evolution operator U of the 
SQRW on the hypercube may be decomposed as U — 
G 2 G X = (1 <g> M)d ...C d , where C a = a x <g> \a)(a\ + 
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1 ® (1 - \a)(a\) is the 0CNOT operator with the tar- 
get qubit being the a-th qubit from the position register, 
and accepting the state \<j>) — \a), with a = l,,..,d. 
The operators C a are mutually commuting, and have the 
common eigenvectors. To find them, we decompose the 
eigenvectors into the product of 2 x d dimensional vectors 
\tp) ® \x)- Applying C a on \ip) |x) we obtain (providing 
that a x \ip) = A|V» 

C a \ip)\x) = <7 s |V)|o)(o|x> + \4>) ® (Ix) - |o)<o|x» 

= M®[(A-l)<a| X >|a> + | X >] 

< (18) 
fl^lx) :A = 1; 

\|V)®[|x>-2<o|x)|o)] :A = -1. 

The case with A = — 1 has to be dealed separately. If 
|x) = \a), then we get C a \ip)\a) = -\ip)\a), and if |x) 
h . b ^ a then C a \ifj)\b) — \i/j)\b). The conclusion: the 
eigensystem of C a is the set of vectors \if>}\x)> where \ip) 
is the eigenvector of a x and x — |a) , a = 1, . . . , d. What 
about the eigensystem of G2? The matrix M can be 
diagonalized, since it is translation invariant. We search 
for the eigenvectors in the form ip^ = J2 a =i d e 2l ™ fea / d e a , 
which yields 

M^ fe = {re 2mka/d + te 2 " k( - aeb)/d )e a 

a=l,d 6=1, d-l 

= (r + t J2 e- 2 " kb/d ) e 2mka / d e a . (19) 

6=1, d-l a=l,d 

where Afe = r — t, if k ^ 0, and A& = r — t + fd 



VII. SQRW IS SUPERSET OF COINED 
QUANTUM RANDOM WALK 

In this section we will discuss the connection between 
the scattering and the coined quantum random walks. 
The SQRW reduces to coined quantum random walk on a 
regular graph (having all vertices with the same degree) , 
and conversely, the SQRW is the generalization of the 
coined quantum random walk on general graphs. 

There is an isomorphism between coined quantum ran- 
dom walk (CQRW) and SQRW on the same Cayley 
graph over Abelian group, Q. We recall that CQRW 
is defined by an unitary operator E on Hilbert space 
He = Hx ® Ha where Hx is spanned by vectors 
\x),x € Gv and Ha is spanned by the generators of Q, 
the basis vectors \a). One step of CQRW is given by 
E = SC where S = £ a T a ® n a and C = 1 <g> M. Here 
T a \x) = \x + a) is the translation, ir a is the projection to 
I a) and M is any unitary operator. The isometry is given 
by 1-1 mapping of basis vectors of both H and He like 
\x(x + a))ji = \x)\a)-n E . The correspondence between 
operators U and E is: U x \yx), where y + a = x, is the 
same like applying the translation S on \x — a)\a) and 



then applying the coin M on \a) s.t. matrix representa- 
tions of U x in natural basis of H x , H x and M in natural 
basis of Ha are the same. 

For regular graphs, we can decompose He into direct 
product of Hex s.t. Hex — span{|x— a) \a) : a is gen. G}. 

The scheme for generalizing coined quantum random 
walk on general graphs was proposed in [l5|. but this 
required an oracle which operate on the set of all edges 
of the graph. Our scheme is based on local operations 
done by multiports, so it is more reasonable and easier 
to implement physically. This was actually proposed in 
|l4j . Algorithms based on coined quantum random walks 
were proposed e.g. in [TJ 

VIII. SEARCHING WITH SQRW 

In this section we will address a question whether it 
is possible to use the SQRW for a database search, or a 
similar task? To answer this question, we need to formu- 
late what a quantum database is and how we can move 
around its entries using the SQRW. 

The database we are searching in is the so called quan- 
tum dictionary. The classical dictionary is a set of pairs 
(key, value). The set of all keys is given by the topology 
of a graph, yielding the adjacency relations among all the 
keys. Random walk (classical) in the dictionary is bound 
to the edges of this graph. 

The searching problem in the dictionary is given as fol- 
lows: given a value, find a key, such that (key, value) is 
in the dictionary. For N keys, this is an O(N) problem. 
To obtain a quantum version of this scheme, we have 
to "quantize" (non-canonically) the problem. Due to the 
fact that the graph is not regular, we cannot factorize the 
complete Hilbert space, but we need to label the states 
in the most general fashion: \x, y), where (x, y) is an 
edge. The searching procedure consists of applying one 
step of the SQRW, and then by querying the database. 
The query corresponds to an application of a unitary op- 
erator (the oracle) |18|. which flips some auxiliary qubit, 
depending on whether the value assigned to the key is the 
one we are searching for. That is, the oracle is the trans- 
formation 0\x,y)\q) h-> \x,y)\q ® f(x,y)}, where f(x,y) 
gives the value 1 if any of the vertices x, y satisfy the 
query, and otherwise. It is clear that the oracle O is 
unitary. Now the searching algorithm is based on the 
sequence of operations (OU) n , where U makes one step 
of the SQRW and O is the oracle query (equivalent to 
the action of the multiports). One such algorithm has 
been presented in Ref. 2|: In our terms it is the SQRW 
on the hypercube, where the multiport assigned to one 
marked key has trivial coefficients r, t (they only change 
the phase), while the other multiports have coefficients 
corresponding to the action of the Grover operator to the 
direction states. In Ref. ,2) it has been shown that the 
marked key can be found in 0(y/~N) steps with probabil- 
ity O(l), where N is the number of the vertices of the 
hypercube. 
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IX. SYMMETRIES OF THE EVOLUTION 
OPERATOR U 

The basic relation Uip = Xtp, where ip = ^l xa lxa\xa) 
yields the following recurrence relation: 



rjx.-a +t^2lxb = ^Ix+a.a 



(20) 



Finding the symmetries of this operator helps us to find 
its eigensystem. We can Fourier transform the states of 
Ti. to another basis, in which the solutions can be easier 
found. The operator U has many symmetries, one of 
them is the translation Tj, : x <— > x + b. The eigenvectors 
of Tb are (for details see Ref. [Tjj) 



|fca)=£(-l) fe *| 



(21) 



with eigenvalues (—l) kb , k E l\. The action of U on 
\ka) is 

U\k~a) = (-l) fe « (r\ka) + J2 \kb}) , (22) 



and in the basis \ka) 
where 



U has the form U = diag ({Vfc}) 



V k = 



fr{-l) kl t(-l) k2 ... t{-l) kd \ 
i(_l)fci r (-l) fc2 t(-l) ks 



(23) 



r (-l) k dJ 



Now we only need to find the eigensystem of this compar- 
atively small matrix. It is obvious, that V k is translation- 
ally sym metric. The eigensystem of (|23f) can be found in 
Ref. [13. 

In what follows we will find another symmetry. Unlike 
the previous case, now we will be changing both the el- 
ements of the position and the direction Hilbert spaces. 
This transformation R will change the vector \x, a) s.t. 
the binary string x is cyclically shifted right by one place 
and a is set to a © 1 modulo d. Since a is unambiguously 



defined by the position in the binary string at which x 
differs from x + a, this transformation is a symmetry. 
This can be viewed of as a rotation about the line seg- 
ment connecting two opposite vertices ... and 1 ... 1. 
We can choose any other two vertices x, y s.t. \x — y\ = d, 
and get a symmetry operator R xy = B\.RB X , where B x 
changes the role of to a; and 1 to 1 + x. More precisely, 
B x \z, a) = \z + x,a) (hence B^ = B). Two transforma- 
tions R x ,R y generally do not commute, but they both 
commute with U. 



X. CONCLUSION 

We have proved that the SQRW is in fact a version of 
the coined quantum random walk. We can use this ob- 
servation to extend the coined quantum random walk to 
the cases of non-regular graphs. While it is in principle 
easy to construct the SQRW on any graph, it is still a 
question whether we also can simulate it efficiently (e.g. 
like in Sec lVI|) . This point is crucial for further devel- 
opment of quantum algorithms based on the SQRW in 
higher dimensions (where the speedup may become no- 
ticeable). The class of algorithms based on SQRW is the 
database searching, using the oracle queries along with 
the "random" steps. We already know at least one such 
algorithm (see Ref. 0) and we know that it is optimal. 
We cannot expect the complexity drop below 0(V~N) for 
N database keys, but the new algorithms may be more 
general in their inputs, and maybe more easy to imple- 
ment. 

We have found the connection between mixing prop- 
erties of the multiport (or the coin), and the distance of 
the respective operator from the unity. It might be in- 
teresting to find an exact function of this distance, which 
yields the measure of mixing for the SQRW. 
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